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Abstract 

The recently derived input-output relations for the radiation field at a disper- 
sive and absorbing four-port device [T. Gruner and D.-G. Welsch, Phys. Rev. 
A 54, 1661 (1996)] are used to derive the unitary transformation that relates 
the output quantum state to the input quantum state, including radiation and 
matter and without placing frequency restrictions. It is shown that for each 
frequency the transformation can be regarded as a well-behaved SU(4) group 
transformation that can be decomposed into a product of U(2) and SU(2) 
group transformations. Each of them may be thought of as being realized by 
a particular lossless four-port device. If for narrow-bandwidth radiation far 
from the medium resonances the absorption matrix of the four-port device 
can be disregarded, the well-known SU(2) group transformation for a lossless 
device is recognized. Explicit formulas for the transformation of Fock-states 
and coherent states are given. 
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I. INTRODUCTION 



Four-port devices such as beam splitters are indispensable to optical investigation, and 
a number of fundamental experiments in quantum optics necessarily require the use of 
them. The quantum theory of dispersionless and nonabsorbing beam splitters has been well 
established |jT|-0. A beam splitter can be realized by a multislab dielectric plate, which is a 
dispersive and absorbing device in general. Even if the effects of dispersion and absorption 
(in a chosen frequency interval) are small, their influence on nonclassical radiation should 
be considered carefully. On the other hand, in practice multislab dielectric configurations 
with strongly varying dispersive and absorptive properties, e.g., near optical band gaps, 
have been of increasing interest, and a description of their action in the quantum domain is 
desired. 

To give a quantum theory of a dispersive and absorbing (linear) four-port device, a 
Kramers-Kronig consistent quantization scheme of the electromagnetic field in dispersive 
and absorbing inhomogeneous media is required PHT^- In particular, quantization of the 
radiation field within the framework of the phenomenological Maxwell theory (with given 
complex permittivity in the frequency domain) can be performed using an expansion of the 
electromagnetic field operators in terms of the Green function of the classical problem and an 
appropriately chosen infinite set of bosonic basic fields This quantization scheme, which 
may be regarded as a generalization of the familiar concepts of mode expansion, applies 
to any inhomogeneous dielectric matter and is consistent with both the Kramers-Kronig 
relations and the canonical (equal-time) field commutation relations in QED |TT|JT^ . 

The formalism has been used in order to derive input-output relations for radiation at 
a dispersive and absorbing (multilayer) dielectric plate and to express the moments and 
correlations of the outgoing fields in terms of those of the incoming fields and the (initial) 
dielectric-matter excitations P, p!0| , p!3| , p!5| . Such a (multilayer) dielectric plate may serve as a 
model for a number of four-port devices, such as beam splitters, mirrors, thin films, interfer- 
ometers, and optical fibers. The results have been used for studying low-order correlations 
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in two-photon interference effects p!3| , [T^JT6| . 



In this paper we extend the input-output relations for the radiation field at a dispersive 
and absorbing four-port device to the complete SU(4) group transformations for radiation 
and matter, and present closed formulas for the transformation of the quantum state as a 
whole. It is worth noting that the theory applies to optical fields at arbitrary frequencies 
and bandwidths. In particular for narrow-bandwidth light in which the frequencies are far 
from medium resonances so that absorption may be disregarded, the well-known results of 
SU(2) symmetry are recognized. In the general case of nonvanishing absorption, for each 
frequency the SU(4) group transformation can be given by a product of eight U(2) and 
SU(2) group transformations, which correspond to an equivalent network of eight lossless 
four-port devices for radiation and matter. 

The paper is organized as follows. In Sec. || the underlying theory is outlined and the 
basic input-output relations are given. The problem of quantum-state transformation is 
studied in Sec. |T1| and closed solutions are presented. To illustrate the theory, exphcit 
transformation rules for Fock states and coherent states are presented. A summary and 
some conclusions are given in Sec. |rV|. 

II. BASIC EQUATIONS 

Let us consider two light beams (of fixed polarization) that propagate along the (positive) 
Xi and X2 axes and impinge on a dispersive and absorbing four-port device that gives rise 
to two outgoing beams propagating along the (positive) yi and 1/2 axes. Following |]T3|, the 



operator of the vector potential in each of the four channels of the device can be given by 







/ duj 




Jo 


A 



^ (1) 



xe'^^(")'""^/^c(2,-,u;) + H.c. 



(j = 1, 2), where 



rij 



H = = PA^) + (2) 
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is the complex refractive index of the adjacent medium on the jth side of the device {A, 
plane area of the beam). In Eq. (|1]), Cj i^Zj , UJ ) stands for the amplitude operators aj{xj,uj) 
and bj{yj, u), respectively, of the incoming and outgoing damped waves at frequency u. The 
input-output relations for the amplitude operators can be derived to be 
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hiVj^^) = J2 Tjjf{uj)af{xf,uj) + J2 Ajj'{^)9j'{^) (3) 
j'=i j'=i 

where it is assumed that the incoming beams enter the device at Xj = Xj and the outgoing 

beams leave the device at yj = yj. The operators gj{uj) play the role of operator noise sources 

and describe device excitations. The 2x2 matrices Tjj/{uj) and Ajji{uj) are the characteristic 

transformation and absorption matrices of the device. Whereas the Tjji matrix describes 

the effects of reflection and transmission, the Ajj/ matrix results from the losses inside 



the device (for Tjji{uj) and Ajjt{uj) of a multilayer dielectric slab, see ||T3|). Finally, the 
commutation rules for the amplitude operators of the incoming waves and the operators of 
device excitations are 

[aj{xj, uo), a],(x^-,, uo')] 

= (5,y(5(cu-a;')e"^''^"^"'"'""^''^', (4) 

[gj{uj),g\,{uj')\ = 6jf6{uj - u'), (5) 

[aj{xj,u),gj,{Lj'y] = (6) 

{xj > Xj). Since the dependence on space of the amplitude operators outside the device is 
governed by quantum Langevin equations, the input-output relations together with the 
commutation relations - (|§) fully determine the action of the device. The commutation 
relations (^) - reveal that the amplitude operators of the incoming waves at the entrance 
plane, aj{uj) = aj{xj,u), and the operators of the device excitations, gj{uj) are independent 
bosonic operators. The amplitude operators of the outgoing waves, bj^u) = bj{yj,u!), do not 
satisfy bosonic commutation relations in general. For given matrices Tjji{uj) and Ajji{(jj), 



their commutation relations can be derived straightforwardly, applying Eq. (|]) and using 
Eqs. ® - (D- 

Let us consider the case when the device is surrounded by vacuum [nj{u!) —>-l]- In this 
case the amplitude operators of the incoming and outgoing waves become independent of 
space and reduce to ordinary bosonic operators. In particular it can be shown that the 
matrix relation 
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(7) 



k=l 



k=l 



is valid, which implies the bosonic commutation relation 

[b,{u):bUuj')] = 5,rS{u-u'). 



(8) 



Note that the relation (|^) reflects the fact that when the device is embedded in vacuum, the 
sum of the probabilities for reflection, transmission, and absorption of a photon is equal to 
one. When the device is embedded in a medium, the matrix relation and the bosonic 
commutation relation (P) are not valid in general. From Eqs. (Q) - (^) it can be seen that a 
unitary transformation 



(9) 



k=l 



[{X ^)jfc = X^J can be introduced such that 



A7^(^)S;.H,A7(a;)6't(a;0 



6jj'6{uj — io') 



(10) 



(Aj >0). Hence the transformed and scaled operators A^ ^ {uj)b'j{uj) are bosonic operators and 
the corresponding (scaled and transformed) transformation and absorption matrices satisfy 
the condition (|^. 

Without loss of generality we can therefore restrict our attention to a bosonic system 
and assume that the matrix relation (|^) is valid. For notational reasons it is convenient to 
introduce the definitions 



a u; 



0,2 (tu) 

5 



(11) 



[UJ 



(12) 



and 



T(uj) 



y ^21 (a;) T22{uj) J 



(13) 



(14) 



Muj) 



(15) 



\^v42l(c^) A22{UJ) j 

The input-output relations for radiation at a general four-port device can then be given in 
the compact form of 



hiuj) = T(uj)ai(uj) + A(uj)g(uj) 



(16) 



with 



T(cj)T+(cj) + A{lu)A+{lu) = I. 



(17) 



III. QUANTUM STATE TRANSFORMATION 

The operator input-output relation (p!6[) enables one to calculate arbitrary correlations of 
the outgoing beams from the correlations of the incoming beams and the device excitations 
13| . To obtain the quantum state of the outgoing beams as a whole, the question arises 



which quantum state transformation corresponds to the operator input-output relation. 
Let us assume that the incoming fields and the device are prepared in a quantum state 
described by the density operator and that for any frequency the input-output relation 
(0) corresponds to the unitary operator transformation 
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b(c.;) = U^k{uj)U, = U-\ (18) 

The effect of the device can equivalently be described by leaving the photonic operators dj{uj) 
unchanged but transforming the input-state density operator to obtain the output-state 
density operator ^out as 

^out = UQinU^. (19) 



A. Lossless device 

Let us first restrict our attention to a field in a sufficiently small frequency interval of 
width Au in which absorption may be disregarded. For this frequency window the four-port 
device can be regarded as being lossless, and Eqs. ([TB|) and (|I7p reduce to 

h{u;) = T(a;)a(a;), (20) 

T(cu)T+(cj) = I. (21) 
Here, the elements of the U(2) matrix T{uj) are usually given by 



( t{uj) r{uj)^ 



giv'H^ (22) 



where t{u)) and r{uj), respectively, correspond to the complex transmittance and reflectance 
of the device at frequency u, 

t{uj) = cos^(cu) e^°("), r = sin^(u;) e'^'-'^K (23) 

When the phase shift (p{uj) can be disregarded, then the U(2) group transformation reduces 
to an SU(2) group transformation. Note that the phase shift can always be included in the 
input operators by replacing sl(uj) with SL(uj)e^'^^^\ so that T(uj) becomes an SU(2) group 
matrix. The unitary exponential operator U in Eqs. (|TB|) and (|I^) can easily be found by 
extending the formalism of lossless beam-splitter transformation |l|JH] to multi-mode fields: 



U = exp 



(24) 



(the superscript T introduces transposition), where the 2x2 Hermitian matrix V(u;) is 
related to the SU(2) matrix T{uj) in Eq. (|20|) as 



(25) 



The operator U can be factored in different ways, e.g.. 



U = exp<^ i / du (f{u)[ai{u)di{u) + a2{u!)d2{^)] 

I JAuj 



X exp 
X exp 
X exp 
X exp 



dujr*{uj) a2(ti;)di(co') 



/ duj lnt{uj) a2{uj)d2{uj) 
Jau 



(26) 



B. Dispersive and absorbing device 



1. Transformation law 



In order to apply the input-output relation (|T6|) [together with Eq. (p]7|)], we first extend 
it to a U(4) group transformation. For this purpose we combine the two-dimensional vectors 
sl{uj) and g{uj) to obtain a four-dimensional input vector 



a (a;) 



^ k(uj) ^ 



(27) 



and supply the two-dimensional vector h{uj) with some other two-dimensional vector h(u;) 
to obtain a four- dimensional output vector 



f3{uj) 



(28) 
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Now we relate the four-dimensional vector f3{uj) to the four-dimensional vector cx{uj) as 

^{lu) = A{uj)a{uj), (29) 



(30) 



where the 4x4 unitary matrix A(co') is chosen such that the input-output relation ([TBD 
between h{uj) and a(ci;) is preserved. The matrix A(ci;) can be expressed in terms of 2 x 2 
matrices as (App. ^ 

/ . ... \ 

A(^)= , (31) 



V 



T(cj) A(c<;) 
-S(uj)C-Huj)T(uj) C(uj)S-Huj)A(uj) 



where 



CM = a/tmt+M 



and 



(32) 



S{uj) = JA{uj)A+{uj) 



are commuting positive Hermitian matrices, and 



C(ujY + S(uj) 



(33) 



(34) 



In Eq. (|31|) the unitary matrix D(c<j) that appears in Eq. ( [All| ) in App. ^ has been omitted, 
since it corresponds to an irrelevant change of the device variables h(w), as it can be seen 
from the second line in the large brackets in Eq. (pij). Note that after separation of phase 
factors e^'^^^^ and e^'^^'^\ respectively, from the matrices T{uj) and A{uj) and inclusion of 
them in the operators sl{uj) and g{uj) the matrix A{uj) can be regarded as an SU(4) matrix. 

The input-output relation (^) can be expressed in terms of a unitary operator transfor- 
mation 



(35) 



where the unitary operator U is given by 



U = exp 



-i / duj(a.'(uj)) ^(uj)a(uj] 
Jo 



(36) 



Here, ^{u) is a 4 x 4 Hermitian matrix which is related to the SU(4) matrix A{uj) by 



exp[-?$(u;)] = A(a;). 



(37) 



Note that for narrow-bandwidth radiation far from medium resonances the u integral in 
Eq. ( P5| ) can be restricted to a small interval in which absorption may be disregarded, A(u;) 



: 0, and hence 



T{uj) 

I 

Y{uj) 





(38) 



(39) 



In this case Eq. (|36D approximately reduces to Eq. (|^) and the SU(2) group transformation 
for a lossless device is recognized. 

Obviously, the first and the second line in the vector equation ( |35| ) correspond to the 
input-output relation (p^. Application of Eq. (O) then yields the output quantum state 



CP) 

Qont, from which the quantum state of the outgoing radiation field, Qo^i, can be derived 



£Jt = Tr(°){^out} = Tr(°){f/^i,m}, 



(40) 



where Tr*-^^ means the trace with respect to the device. The input density operator is 
an operator functional of cx.{uj) and a.^{uj), 



and hence the transformed density operator ^out can be given by 



(41) 



&ut = ftn[f/«(^)^7^^7At(^)m]. 



(42) 
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Recalling Eqs. p9|) and (pS]), we see that 

Uct{uj)U^ = 



(43) 



(44) 



Combining Eqs. (^OD - (|4^), we derive 

^iul = Tr(°){ft„[A+(^)a(^), A^(.;)«t(^)]} 



(45) 



5. Relation to U(2) and SU(2) group transformations 

As shown in App. |B|, the U(4) group transformation defined by the matrix A given 
in Eq. (|A11|) is equivalent to five U(2) group transformations. That is to say, for chosen 
frequency the action of an absorbing four-port device formally corresponds to the combined 
action of five lossless four-port devices in general (for the factorization of an U(N) matrix 
into U(2) matrices, see also jl^)- When the irrelevant matrix D(c(j) in Eq. ([All|) is set equal 
to the unit matrix I, then Eq. ([All|) reduces to Eq. ( pT)) and the action of the absorbing 
device corresponds, for chosen frequency, reduces to the combined action of eight lossless 
devices in general. In this case the unitary operator (EB) can be factored into a product of 



unitary operators of the type given in Eq. (|2lD for a lossless device. 



IJ\M] q] = exp 



doj (q^(cu))'^W(cu)q(cj) 



(46) 



Here, W(a;) is a 2 x 2 Hermitian matrix that is related to a U(2) group transformation 
matrix M(a;) as 



exp[-zW(cj)] = M(tu) 



(47) 



and q(co') is a vector whose two components are bosonic operators. Note that for narrow- 
bandwidth radiation far from medium resonances Eq. ( ^61) [together with Eq. (^T])] corre- 
sponds to Eq. (^ [together with Eq. (^], with M(cj) = T(cj), W(cj) = V(cj), and q(cu) = 
k{uj). As shown in App. 0, the unitary operator given in Eq. (|36|) , 
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U = U[A; a] = exp 



Jo 



(48) 



can be decomposed into a product of operators [/[M; q] as follows: 

U[A; a] = U[C+iS- {ik+g)/V2] 

xU[C-tS; (a+?g)/V2] U[S-^A; g] U[C-^T;§i] 



(49) 



[cf. Eqs. (^), ( PTSD , and (iH)], and decomposition of f/[C-zS; (a+ig)/v^] and U[C + 
iS; (ia+g)/-\/2] eventually yields 



f7[A; a] 



X U[C+tS;g]U[C-tS;Si] 

x?7[P;d2]f/[P;di] 

x?7[S-iA;g]f/[C-iT;a], 



(50) 



where 



aj{uj) 



(51) 



(j = l,2) and 



(52) 



[cf. Eqs. ( P26| ) and (P27|) ]. It should be pointed out that when A{u)) is an SU(4) group 
transformation, then the matrices P, S^^(co')A(ci;), and C^^{uj)T{uj) correspond to SU(2) 
group transformations. The matrices C(ti;) +28(^7) and C(ci;) — iS{uj) correspond to U(2) 
group transformations in general, i.e., SU(2) transformations and additional phase shifts. 
Needless to say that each of the operators f/[M; q] on the right-hand side in Eq. ( |50D can 
be further factored, e.g., according to Eq. (p6|). 
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3. Discretization 



In quantum optics radiation fields are frequently described in terms of discrete modes. 
Here we restrict attention to (quasi-)monocliromatic discrete modes. For this purpose we 
subdivide the frequency axis into sufficiently small intervals with midfrequencies Um and 
define the bosonic input operators 

otm = —r^= [ duj a{uj) (53) 



and the bosonic output operators accordingly. The operator input-output relation (p9D 
then reads as 

Pm = ^mOCm, (54) 

[A^ = A{um)], which can be rewritten as 

= U^dCmtj = Ul^OLmUm , (55) 



where [in place of (^)] 



U = JlUm. (56) 



with 



Urn = eKv[-K^lnf^rnOLrn\, (57) 

and according to Eq. (|37D , the 4x4 Hermitian matrix is related to the SU(4) matrix 

exp[-z$^] = A^. (58) 

The input density operator is now an operator function of ct^ and a)^, and according 
to Eq. (45), the density operator of the outgoing radiation field can be given by 

£1 = Tr(°){ft,[A+a^, Xl&l]} . (59) 
13 



In close analogy to Eqs. (|49|) and (0), each SU(4)-group transformation operator Um can be 
decomposed into a product of U(2)- and SU(2)-group transformation operators Um[^m] Qm], 

Um[-^m] f^m] = Um[Cm-\-iSm] (^^r?i + §m)/V^] 

xf/[S„iA„;gJf/[C„iT^;aJ 
xf/[P;d„2]f>m[P;d„i] 



(60) 



Here f/m[M„; q^] is given by 



Um[Mrn; Qm] = CXp -«(q|„)^W„4 



(61) 



where 



exp[-zW™] = M, 



(62) 



[cf. Eqs. ( ^61) and (|i7|)]. Recalling the definition of discrete operators, Eq. (|53|), application 
of Eq. ( P^ to f/m[Mm; Qm] for further factorization is straightforward. 



4- Transformation of Fock- states and coherent states 

To illustrate the theory, let us consider the transformation of Fock states and coherent 
states as fundamental basis states for quantum state representation. For the sake of trans- 
parency, we will restrict attention to single-mode states at some chosen frequency, so that 
the subscript m can be omitted. The results can easily be extended to multimode fields by 
taking the direct product of single-mode states. Let be 

0m = \^m){^m\, (63) 
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lipin) = \ni; na, ng; n^) = J| —7=^ |0), (64) 



^=1 Vn^! 

the density operator of the system in the case when rii and n2 photons impinge on the device 
that is excited in Fock states with and quanta. From Eq. ( pQ] ) we then obtain 

£l = Tr(°){|^out)(^out|}, (65) 

with 

1^0,0 = n^ IE v«ii io)- (66) 

We use the decomposition 



EA>i| = E (67) 



where the (non-negative) integers ku/j. satisfy the condition 

4 



and arrive at 



E ^^/^ = '^'^ ' (68) 



l^out) = E C'fci,fc2,fc3,fc4l^i; ^2; A;3; /C4) (69) 



where 



ki,k2,k3,k4 



= n\/^ nv^ E n (70) 

the ky^ satisfying the conditions 

4 4 

^i^/i = "^v 1 = • (71) 

Using Eqs. (|65|), (|6^), and ([70|), the quantum state of the outgoing radiation-field modes can 
easily be obtained, 

£1= E E ^fci,fe2,fc'i,fc^, l^i;^2)(A:'i;fc2l, (72) 
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(73) 



The density operator of the outgoing field, Qo^^, can be represented in another way which 
more clearly shows the infiuence on the outgoing field state of the device. Let us define 
linear combinations 



(74) 



and 



yl = Y.A2+ru9\ (75) 
of the photonic and device operators, respectively. Making in Eq. (|66D the insertion 

4 

K>^^1 = ^i + yl, (76) 

11=1 

Eq. (|65D then reads as 



(77) 



where 



n 



1 fn,. 



1 V^'^, 



n 



' n 



1 ^ nj 



X 



(78) 



Vi/=1 



and |0*^'^^) (|0*^'^^)) is the ground state of the field (device). The device vacuum expectation 
value in Eq. (|78|) can be calculated by moving the operators y^, from left to right and 
employing the commutation relations between tji, and yj^. Then we find, that the coefficients 
y{p^}{q^} = Y{Pt,}{qv}{.Zun) are functions of the matrix elements Z^^ of the matrix 



^I-T+T -T+A ^ 



V 



(79) 



AT I-A+A 

Finally let us consider the transformation of coherent states. 
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Wm) = |7i;72,73;74) 



Yl exp{-f^al - -fla,y^ |0). 



(80) 



From Eq. (^) we again obtain Eq. ([65|), where l^/'out) is the coherent state 



with 



l^out) = |Ai; A2, A3; A4) 



n exp(A^dJ. - A;;a^) |0), 



(81) 



A^ = X! 



(82) 



From Eqs. (^) and (^) it follows that the outgoing modes are prepared in coherent states, 



^iut — |Ai; A2)(Ai; A2I. 



(83) 



Note that when the device is excited in a coherent state, then the coherent amplitudes Ai 
and A2 of the outgoing modes are not only determined by the characteristic transformation 
matrix T but also by the absorption matrix A and the coherent-state amplitudes of the 
device, as it can be seen from Eq. (|82D , 

2 



Ai = E i'^iJ^j + ^ii7i+2) 



(84) 



(^ = 1,2). 



IV. SUMMARY 

We have developed a quantum theory of the action of a dispersive and absorbing optical 
four-port device, such as a beam splitter. In particular we have presented formulas for calcu- 
lating the complete quantum state of the outgoing fields from the input quantum state of the 
incoming fields and the device excitations, without any frequency restriction. The theory 
is a natural extension of the standard theory of lossless beam splitters. According to the 
underlying quantization scheme for radiation in inhomogeneous Kramers-Kronig media, the 
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device is described in terms of a frequency- dependent transformation matrix that includes 
transmission and reflection and a frequency-dependent absorption matrix. 

For each frequency the action of the device has been described in terms of a U(4) group 
transformation of incoming field operators and device operators. Each U(4) group transfor- 
mation can be realized in a natural way by the combined action of eight lossless four-port 
devices. However, each U(4) matrix is only determined up to a U(2) matrix. This matrix 
can be chosen such that the U(4) group transformation is equivalent to five U(2) group 
transformations. That is to say, for chosen frequency the action of an absorbing four-port 
device formally corresponds to the combined action of five lossless four-port devices. 

The quantum state of the outgoing radiation can be expected to sensitively depend on 
the quantum state the device is prepared in when the incoming fields impinge on the de- 
vice. In combination with conditional measurement this offers novel possibilities of quantum 
state manipulation. In particular, the theory enables one to study the effect of resonance 
frequencies on quantum state transformation. 
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APPENDIX A: DERIVATION OF THE U(4) GROUP MATRIX 



Let us write the sought U(4) matrix A(co') as 

/ 

A(a;) = 



T(cj) A(cj) 



(Al) 



where T{uj) and A{uj) are defined in Eqs. ( |14|) and (|T5D and satisfy the relation ([T7|). The 
2x2 matrices F{uj) and G(ci;) are to be determined such that A{uj) is unitary, i.e.. 



FMF+M + GMG+M = I, 



(A2) 



F(w)T+(cj) + G(cu)A+(cu) = . 



(A3) 



From Eq. (^^) we find that 



We substitute in Eq. (|A2|) for F{u) the result of Eq. (|A4|) and derive 



(A4) 



G(cu) I + A+(lu) T{u)T+{uj) A(cu) G+(a;)=I 



(A5) 



and hence 



I + A+(cu) THT+(cj) AH = G+(cu)G(cu 



-1 



(A6) 



Recalhng Eq. (|T7D, from Eq. ( |X6| ) we find that 



G+HG(cj) = I - A+(cu)AH . 



(A7) 



A particular solution of Eq. (|A7|) is 



G{uj) = C(a;)S-^(a;)A(a;), (A8) 

where C(ci;) and S(uj) are defined in Eqs. (|32D and (^31), respectively. Obviously, the general 
solution reads as 



(A9) 



where D is an arbitrary unitary 2x2 matrix. From Eq. ( [A^ ) it then follows that F{uj) is 
given by 



F{uj) = -B{u)S{uj)C-\uj)T{lu). 



(AlO) 



T(a;) 



(All) 



Combining Eqs. (|A1|) , ([A9|), and (|A1CI|) , we obtain 
/ 

A(.;) = 

^ -D(o;)S(o;)C-nc^)T(cj) F>{uj)C{uj)S-\uj)A{uj) 

which reveals that for given matrices T{uj) and A{uj) the U(4) matrix A{uj) is only deter- 
mined up to a U(2) matrix F){uj). 



J 
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APPENDIX B: FACTORIZATION OF THE U(4) GROUP TRANSFORMATION 



The U(4) matrix A(co') in Eq. ( All ) can be rewritten as a product of three U(4) matrices 
as follows: 



A{u) = Asiuj) A2(^) Ai(a;), 



(Bl) 



where 



Ai(^) 











DMS-^MAM 



(B2) 



A2(^) 



^ D(cu)C(cj)D+(cj) DHS(cu)D+(u;) ^ 



-D(cj)SHD+(cj) D(cj)C(cu)D+(cj) 



(B3) 



( DtM 0^ 



V 



I 



(B4) 



When we choose the matrix D(ti;) such that the matrices D(u;)C(co')D+(co') and 
D(co')S(ti;)D~'"(ti;) become diagonal matrices [note that C(a;) and S(u;) defined in Eqs. (|32D 
and (0), respectively, can be diagonalized by the same unitary matrix], then the U(4) group 
transformation corresponds to five U(2) group transformations. 

Let D(a;) be the unit matrix, D(a;) = 1. In this case Eq. ( |BI] ) reduces to 



A(lu) = A2(uj) AUco) 



(B5) 



where now 



^ C-Ulj)T(uj) 



Aiicu) 







^ s-^HAH^ 



(B6) 



and 



20 



V 



(B7) 



The matrix S.2[iS) can be given by the unitary transform of a quasi- diagonal matrix ^^^^^ 



A2H = T+A2(cj)T 



(B8) 



where 



C(cj) -zS(cu) 







c(cj) + iSH 



(B9) 



and 



(BIO) 



Combining Eqs. (p5| ) and (p8|), we obtain 



A(cu) = T+A'2(cj)TAi(cu) 



(Bll) 



which corresponds to a decomposition of the U(4) group transformation into eight U(2) 
group transformations. 

Using Eqs. (^9|) and (|B5D and recalhng Eqs. (|35D - (|3^) , we may write 



A2(c<j) Ai(c<j) diiuj) 
A2(a;) [/|d(u;)f/i 
f/| AaMAHf/i 



(B12) 



with 



f/ = ?7[A; a] = [/[Aa; a]f/[Ai; a] 



Here, t/j = f/[Ai; a] (i = 1, 2) is given by Eq. (|36|) , with $i(ct;) in place of and 



(B13) 
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exp[-i$i(cj)] = Ai(cj). 



(B14) 



From the quasi- diagonal structure of Ai(u;), Eq. (|B6|) , it then follows that 

^Wi(cj) 0^ ( ^ 




\ 



J 



(B15) 



\0 W2H J 



where 



(B16) 



and 



(B17) 



Thus, U[Ai; a] can be expressed in terms of two unitary operators of the type given in 
Eq. (IID [together with Eq. 

f/[Ai; a] = U[S-'T- g] U[C-'T; a]. (B18) 

To decompose U[A2; ct], we note that Eq. ( [B8D implies that 

?7[A2;a] = f/[A'2; T(cu)a], 



(B19) 



where 



exp[-z$2(^)] = ^2(^) 



(B20) 



and 



(B21) 



The quasi-diagonal structure of A'2{uj), Eq. (P^D, enables us to write 





0^ 


+ 




^ 


V 








-W3(l^) ^ 



(B22) 
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where 



expl-iWsiuj)] = CH -zS( 



(B23) 



so that U[A2] ct] can also be expressed in terms of two unitary operators of the type given 
in Eq. [together with Eq. (g^)], 



U[\2; a] = U[C+tS; itk+g)/V2] 



(B24) 



Recalhng the definitions of dj^u), Eq. (0), and P, Eq. (^), it is seen that (a + ig)/-\/2 and 
(isL + g)/ can be given by 



1 ( o-jiuj) +igj{io) ^ 
72 



Pdj-(cu) 



= [/t[P;d,]d,Hf/[P;d,] 



(j = 1, 2), and hence 



U[C-tS; {k+zg)/V2] = f/t[P; d^] f/t[P; d^] 

X [/[C-iS; a] J7[P; da] f/[P; di] 



(B25) 



(B26) 



U[C + iS; {ik+g)/V2] = f/t[P; di] U^[P; ds] 

xf/[C-2S;g]?7[P;d2]f/[P;di]. 



(B27) 
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